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Abstract 

Dynamically triangulated surface model is found to undergo a first-order crumpling 
transition between the smooth phase and a crumpled phase. Triangulated spheres 
are divided into compartmentalized domains, whose boundary bonds remain un- 
flipped in the Monte Carlo simulations. The vertices can freely diffuse only inside 
the compartments, and they are prohibited from the free diffusion over the surface 
due to the domain boundaries. The violation of reparametrization invariance of the 
Hamiltonian is a reason for the existence of the long-range order in the fluid surface 
model. 
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1 Introduction 



The crumpling transition of membranes is an interesting topic in the softmat- 
ter physics as well as in the biological physics [1,2,3]. A well-known model for 
such transition is the surface model of Helfrich, Polyakov, and Kleinert [4,5,6]. 
A considerable number of theoretical and numerical studies have been devoted 
to reveal the phase structure of the model [7,8,9,10,11,12,13,14,15,16,17]. Re- 
cently, it was shown by Monte Carlo (MC) simulations that the model un- 
dergoes a first-order transition on spherical and fixed connectivity surfaces 
[11,12], and the transition is universal [13]. The vertices can move only locally 
on the surface because of the fixed connectivity nature in those surface models. 
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However, the crumpling transition is not yet clearly understood in biological 
membranes. If we consider the possibility of the transition in the cell mem- 
branes, we should take account of the fluid nature such as the lateral diffusion 
of lipids. 

Conventionally, the free diffusion of lipids has been realized by the dynami- 
cal triangulation technique in the surface model[14,15, 16,17]. The diffusion of 
lipids corresponds to the reparametrization invariance of the model, because 
the diffusion has no cost in energy on fluid surfaces. 

In the cell membranes, however, the free diffusion of lipids is suffered from 
heterogeneous structures. Lipid molecules of the cell membranes are known 
to undergo the so-called hop diffusion over the surface, which was recently 
found experimentally [18]. The free diffusion of lipids is prohibited due to 
the cytoskeleton. The diffusion rate is, therefore, 10-100 times lower than 
that of artificial membranes, which is usually homogeneous and has no such 
domain structure. Moreover, some artificial membranes are considered to have 
skeletons, because they are partly polymerized [19]. 

Motivated by this fact observed in the cell membranes, we study in this Letter 
a dynamically triangulated surface model with compartmentalized domains 
whose boundaries are consisting of triangle edges (or bonds) that are not to be 
flipped. The diffusion is constrained so that vertices can diffuse only inside each 
domain, and hence the vertices never jump across from one domain to the other 
domains. Nevertheless, we consider that such constrained lateral diffusion can 
simulate the hop diffusion in the cell membranes as a first approximation. 



2 Model 



By dividing every edge of the icosahedron into L-pieces of the uniform length, 
we have a triangulated surface of size N = 10L 2 + 2 (= the total number of 
vertices). The starting configurations are thus characterized by N 5 = 12 and 
N 6 = N — 12, where N q is the total number of vertices with the co-ordination 
number q. 

The compartmentalized structure is built on the surface by keeping the bound- 
ary bonds unflipped in the MC simulations with dynamical triangulation. The 
boundary of the compartment is constructed from a sequence of bonds that 
remain unflipped. The total number Nc of the compartments depends on the 
surface size N. We fix n the total number of vertices inside a compartment to 
the following three values: 

ra = 21, 36, 66 (# of vertices in a compartment). (1) 
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As a consequence, Nc is increased with the increasing N. The reason why we 
fix n is that the size of compartment is considered to be finite, and then it is 
expected that total number of lipids in the compartment also remains finite in 
the cell membranes. We must emphasize that the finiteness of n is physically 
meaningful rather than the value of n, because we do not always have one to 
one correspondence between the vertices and the lipid molecules. 

Figures l(a),(b) show surfaces of (N, n) = (2562, 21) and (N, n) = (15212, 66) 
for the starting configurations of MC simulations. Thick lines denote the com- 
partment boundary consisting of the bonds that are not to be flipped. Ver- 
tices on the boundary of compartments can locally fluctuate, and they are 
prohibited from the diffusion. The other vertices freely diffuse only inside the 
compartment. 




(a) (iV,n) = (2562, 21) (b) (N, n) = (15212, 66) 

Fig. 1. Starting configuration of surfaces of (a) (N,n) = (2562,21) and (b) 
(N,n) = (15212,66), where n is the total number of vertices inside one compart- 
ment. Thick lines denote the compartment boundary consisting of the bonds not to 
be flipped. 

We note that the tethered surface model is obtained from the compartmen- 
talized model in the limit n — > 1, where the vertices are prohibited from the 
free diffusion. On the contrary, we obtain the fluid surface model in the limit 
of n^N, where all the vertices freely diffuse over the surface. 

The surface model is defined by the partition function 

Z = J2 I n^expf-S^T)], (2) 

T i=l 

S(X,T) = S 1 + bS 2 , 

where b is the bending rigidity, /' denotes that the center of the surface is 
fixed in the integration. S(X, T) denotes that the Hamiltonian S depends on 
the position variables X of the vertices and the triangulation T. Y/r denotes 
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the sum over all possible triangulations T, which keeps the compartments 
unflipped. The Gaussian term Si and the bending energy term S 2 are defined 
by 

S 1 = J2(X i -X j ) 2 , S 2 = £(1-^.11,), (3) 

(ij) (ij) 

where J2(ij) in 5*1 is the sum over bonds (ij) connecting the vertices i and j, 
and J2(ij) m S2 is also the sum over bonds (ij), which are edges of the triangles 
i and j. rii in Eq. (3) is the unit normal vector of the triangle i. 



3 Monte Carlo technique 

The vertices X are shifted so that X' = X + 5X , where 5X is randomly cho- 
sen in a small sphere. The new position X' is accepted with the probability 
Min[l,exp(-AS)], where AS = S (new) - S (old) . 

The summation over T in Z of Eq.(2) is performed by the standard technique 
for bond flip. The bonds are labeled with sequential numbers. The total num- 
ber of bonds Nb is given by Nb = 3N — 6, which includes the bonds making 
the boundary of compartments. Firstly, the odd-numbered bonds are sequen- 
tially chosen to be flipped, and secondly the remaining even-numbered bonds 
are chosen. The flip is accepted with the probability Min[l, exp(— AS*)]. In 
this procedure, the compartment boundary remains unflipped. iV-updates for 
X and Nb/2- updates for T are consecutively performed and make one MCS 
(Monte Carlo Sweep). The radius of the small sphere for 5X is chosen so that 
the rate of acceptance for X is about 50%. We introduce the lower bound 
1 x 10 -8 for the area of triangles. No lower bound is imposed on the bond 
length. We use surfaces of size iV = 2562, iV = 5762, and iV = 10242 for n = 21, 
and N = 1002, N = 4002, and N = 9002 for n = 36, and N = 1692, N = 6762, 
and N = 15212 for n = 66. 



4 Results 

The mean square size X 2 is defined by 

x2 = ^T,{x l -x)\ X = ^X U (4) 
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Fig. 2. The mean square size X 2 vs. b obtained on the surfaces of (a) n = 21, (b) 
n = 36, and (c) n = 66. Solid lines were obtained by multihistogram reweighting 
technique. 

where X is the center of the surface. Figures 2(a), (b), and (c) show X 2 ob- 
tained on the surface of n — 21, n — 36, and n = 66, respectively. Solid lines 
connecting data were obtained by multihistogram reweighting technique [20]. 
Discontinuous change of X 2 can be seen in each case of n when the size N in- 
creases, and indicates a first-order transition. We also find that the transition 
point b n moves left on the b axis as n decreases. It is expected in the limit of 
n — > 1 that b n reduces to the value, which corresponds to the transition point 
of the tethered surface model [12]. On the contrary, b n moves right on the 
b-axis as n increases. We expect that b n disappears in the limit of n — > N at 
sufficiently large N, which will be confirmed later in this Letter. 
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Fig. 3. The bending energy S2/NB vs. b obtained on the surface of (a) n = 21, (b) 
n = 36, and (c) n = 66. Nb(=3N — 6) is the total number of bonds. 

The bending energy S 2 /N B defined in Eq.(3) is shown in Figs. 3(a), (b), and (c), 
which correspond to n — 21, n — 36, and n = 66, respectively. Discontinuous 
change of 62 /iV^ can also be seen at b where X 2 discontinuously changes, 
although the discontinuity of S2/NB is not sufficiently clear in the figures. We 
find also that the value of S 2 /N B in the smooth phase (or in the crumpled 
phase) increases as n decreases at the transition point. The value of S 2 /N B 
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at the transition point becomes identical to the one of the tethered surface 
model in the limit of n— > 1 [12]. On the contrary, the discontinuous change of 
S2/NB is expected to disappear at sufficiently large n, because the first-order 
transition at n = 21 ~ 66 disappear when n^N at sufficiently large N. 




Fig. 4. The variation of the gap of S2/NB against n, which are obtained on the 
(N,n) = (10242,21), (N,n) = (9002,36), (N, n) = (15212,66) surfaces. The solid 
circles (•) and the dashed lines denote the value of S2/NB in the smooth (the 
crumpled) phase at the first-order transition point. The symbol ((g)) is unknown 
critical value of n denoted by n c , where the discontinuous transition is expected to 
disappear. The position of n c is also unknown on the 6-axis. 

This implies that there exists a finite n c , where the first-order transition turns 
to a continuous or a higher-order one. Figure 4 shows the variation of the gap 
of S2/NB against n. The gap of S2/NB at n = 21 reduces as n increases and 
eventually goes to zero at n = n c , which will be clarified in the future. 

We note that the maximum co-ordination number g max is as follows: g max = 38 
on the (N, n) = 10242 surface at 6= 1.38, g max = 42 on the (N, n) =9002 surface 
at 6 = 1.42, and g max = 40 on the (N, n) — 15212 surface at 6=1.51. These were 
obtained throughout the MC simulations. 

In order to show the discontinuity in S2/NB more clearly, we plot in Figs. 5(a)- 
5(f) the distribution (or histogram) h(S2) of S2/NB and the corresponding 
variation of S^/iVe. They were obtained on the surfaces of (N, n) = (5762, 21), 
(4002,36), and (6762,66). The discontinuity of S2/NB can be seen in the 
histogram on the iV>5762 surfaces. Because of the size effect, the transition 
appears to be continuous on the iV = 4002 surface. The double peak at (N, n) = 
(5762,21) is more clear than that at (N,n) = (6762,66), because the gap of 
S2/NB reduces as n increases. We should note that a double peak structure is 
very hard to see in /^(S^) on the iV > 9002 surfaces. When the configuration 
is once trapped in the smooth (the crumpled) state it hardly changes to the 
crumpled (the smooth) state at the transition point on such large surfaces. 
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Fig. 5. The histogram h(S<i) and the corresponding variation of S2/NB against MCS 
obtained on the surfaces of (a),(b) (N,n) = (5762,21), (c),(d) (N,n) = (4002,36), 
and (e),(f) (N, n) = (6762, 66). 

This problem may be resolved with more sophisticated MC techniques than 
the current one in this Letter. 
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Fig. 6. (a) X 2 vs. b, (b) S2/NB vs. b, and (c) Cs 2 vs. b obtained on the fluid surface 
model without the compartmentalized structure. The symbol n = N drawn on the 
figures denotes that the surfaces are those without the compartment. 

The phase transition is expected to disappear from the fluid surface model 
defined on the surfaces that have no compartment. In order to show this, we 
performed MC simulations on the surfaces without the compartments up to 
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the size N = 4842. Figure 6(a) shows X 2 against b. We can see no abrupt 
growing of X 2 in the figure. The bending energy ^/-/Vg shown in Fig. 6(b), 
where the variation of S2 against b seems almost independent of N. The specific 
heat Cs 2 , which is defined by Cs 2 = j^( (S2 — (S 2 )) ), is expected to reflect the 
phase transition. However, we can see no anomalous behavior in Cs 2 shown in 
Fig. 6(c); there can be seen no peak in Cs 2 - Thus we confirmed that the phase 
transition disappears from the model if n — > iV at sufficiently large N. 




(a) iV= 15212, n = 66 (b) the surface section 

Fig. 7. Snapshot of the surface of (N, n) = (15212, 66) obtained in the smooth phase 
at (a) b = 1.53, which is close to the transition point, and (b) the corresponding 
surface section. 

Finally, we show in Fig. 7(a) a snapshots of the (N,n) = (15212,66) surface, 
which is obtained in the smooth phase close to the transition point. Figure 
7(b) shows the surface section. Thus we confirmed that the smooth phase can 
be seen close to the transition point at finite b in contrast to the case n— >iV, 
where it can be seen only at b— >oc. 



5 Summary and Conclusion 

To the summary, we have studied a compartmentalized surface model for the 
cell membranes and found that the model undergoes a first-order transition 
between the smooth phase and the crumpled phase. The model is classified as 
a fluid surface model, although the long-range order and the phase transition 
can be seen at finite b. The compartmentalized structure is an origin of a 
broken symmetry for the reparametrizaition invariance, which is considered to 
be a reason for the existence of the phase transition. Consequently, the result is 
not in contradiction with the standard argument for the non-existence of long- 
range order in fluid membranes. Moreover, the critical point of the transition 
is strongly expected at finite n c (> 66). The number n c and the corresponding 
critical bending rigidity b c remained to be clarified. However, it is possible that 



8 



the compartmentalized fluid surface model undergoes a continuous transition, 
because the strength of the transition is expected to be continuously softened 
with the increasing n. We observed that the phase transition remains first- 
order at n = 66 and disappear at n^N in the limit of N-^oo. Although the 
demonstration of the continuous transition is computationally a hard task, it 
will be performed and the results including those in this Letter will be reported 
elsewhere. 

This work is supported in part by a Grant-in- Aid for Scientific Research, No. 
15560160. 
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